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Abstract 

An analysis of electroweak corrections to the anomalous dipole moments of lepton from some 
special two- loop diagrams where a closed neutralino/chargino loop is inserted into relevant one-loop 
diagrams of the standard model is presented in the split supersymmetry scenarios. Considering 
the translational invariance of the inner loop momenta and the electromagnetic gauge invariance, 
we get all dimension 6 operators and their coefficients. After applying equations of motion to the 
external leptons, we obtain the anomalous dipole moments of lepton. The numerical results imply 
that there is parameter space where the contribution to the muon anomalous magnetic dipole 
moment from this sector is perhaps significant, and the contribution to the electron electric dipole 
moment from this sector is sizable enough to be observed in next generation experiments. 
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I. INTRODUCTION 



At both aspects of experiment and theory, the magnetic dipole moment (MDM) of lep- 
ton draws the great attention of physicists because of its obvious importance. The anomalous 
dipole moments of lepton not only can be used for testing loop effect in the standard model 
(SM), but also provide a potential window to detect new physics beyond the SM. The current 
experimental result of the muon MDM is [l| 

a e * p = 11 659 208 ± 6 x 1(T 10 . (1) 

From the theoretical point of view, contributions to the muon MDM are generally divided 
into three sectors 2j: QED loops, hadronic contributions and electroweak corrections. With 
the hadronic contributions which are derived from the most recent e + e~ data, we can get 
the following SM predictions jj], ^,0]: 

a SM = 11 659 180.9 ± 8.0 x 10~ 10 , 
a SM = 11 659 175.6 ± 7.5 x 10~ 10 , 



a SM = n g59 17QA ± 9 3 x 10 -io _ / 2 ) 



The deviations between the above theoretical predictions and the experimental data are all 
approximately within error range of ~ 2a. Although this ~ 2a deviation cannot be regarded 
as strong evidence for new physics, along with the experimental measurement precision and 
theoretical prediction accuracy being constantly improved, this deviation may become more 
significant in near future. 

In fact, the current experimental precision (6 x 10~ 10 ) already puts very restrictive bounds 
on new physics scenarios. In the SM, the electroweak one- and two-loop contributions 
amount to 19.5 x 10~ 10 and —4.4 x 10~ 10 [(J respectively. Comparing with the standard 
electroweak corrections, the electroweak corrections from new physics are generally sup- 
pressed by A^ /AjL, where A EW denotes the electroweak energy scale and A NP denotes the 
energy scale of new physics. 

Supersymmetry (SUSY) has been considered a most prospective candidate for new physics 
beyond the SM. Nevertheless, the softly broken SUSY at electroweak scale induces many un- 
wanted phenomenological problems, such as new sources of flavor changing neutral currents, 



new CP violating phases etc. In order to solve those problems, the authors of literature [?| 
have recently proposed a split scenario. In this split scenario, SUSY is broken at a high 
energy scale that could be even near the scale of grand unification theory (GUT). A di- 
rect result of this assumption is that the scalar superpartners of SM fermions are all super 
heavy. On the other hand, charginos and neutralinos acquire masses around electroweak 
scale to TeV or so because of R-symmetry and PQ symmetry. Within this framework, 
heavy sfermions suppress the one-loop supersymmetric corrections to the processes of flavor 
changing neutral currents and lepton anomalous dipole moments at a negligible level. The 
leading contributions of new physics to theoretical predictions only arise from the two-loop 
diagrams in which a closed neutralino/chargino loop is inserted into relevant one loop SM 
diagrams, and the corresponding theoretical corrections are sizable enough to be well within 
the sensitivity of the next generation of experiments 8|. 

Actually, the two-loop electroweak corrections to the anomalous dipole moments of lepton 
are discussed extensively in literature. Utilizing the heavy mass expansion approximation 
(HME) together with the corresponding projection operator method, Ref. 9| has evaluated 
the two-loop standard electroweak corrections to the muon MDM. Within the framework 



of CP conservation, the authors of Ref. 
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11| present the supersymmetric corrections 



from some special two-loop diagrams where a close chargino (neutralino) or scalar fermion 



loop is inserted into those two-Higgs-doublet one-loop diagrams. Ref. [12] discusses the 
contributions to the muon MDM from the effective vertices H ± W Zf, ~f, h (H Q )'~f'-f which are 
induced by the scalar quarks of the third generation in the minimal supersymmetric extension 



of SM. In the split scenario, the supersymmetric contributions to 
(EDM) of lepton have been already presented in Ref. 
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;he electric dipole moment 



141 ] . Under the assumption 



l m i,2|, l/fffl ^ m w! 1 the authors derive 5 CP-odd dimension 6 operators involving gauge 
fields and Higgs after they integrate out charginos and neutralinos at one-loop level. Inserting 
the effective couplings from those CP-odd dimension-6 operators into those relevant SM one- 
loop diagrams, they then get the lepton EDMs. 



1 m 2 , nil denote the masses of SU{2) x U(l) gauginos in the soft breaking terms, and \i H denotes the 
/i-parameter in the superpotential, respectively. 
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FIG. 1: The two-loop self energy diagrams which lead to the lepton MDMs and EDMs in split 
SUSY, the corresponding triangle diagrams are obtained by attaching a photon in all possible ways 
to the internal particles. In concrete calculation, the contributions from those mirror diagrams 
should be included also. 



In this paper, we apply the effective Lagrangian method to get the anomalous dipole 
moments of lepton. The effective Lagrangian method has been adopted to calculate the 
twc-.oop sopersynunetric correction^ to the branching ratio of 6 - SJ Q, neutron EDM 
[1 61 ] and lepton MDMs and EDMs [17]. In concrete calculation, we assume that all external 
leptons as well as photon are off-shell, then expand the amplitude of corresponding triangle 
diagrams according to the external momenta of leptons and photon. Using loop momentum 
translational invariance, we formulate the sum of amplitude from those triangle diagrams 
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which correspond to the corresponding self-energy in the form which explicitly satisfies the 
Ward identity required by the QED gauge symmetry. Then we can get all dimension 6 oper- 
ators together with their coefficients. After the equations of motion are applied to external 
leptons, higher dimensional operators, such as dimension 8 operators, also contribute to the 
muon MDM and the electron EDM in principle. However, the contributions of dimension 8 
operators contain an additional suppression factor m 2 /A 2 p comparing with that of dimen- 
sion 6 operators, where mi is the mass of lepton. Setting A NP ~ lOOGeV, one finds that this 
suppression factor is about 10~ 6 for muon, and 10~ 10 for electron separatively. Under current 
experimental precision, it implies that the contributions of all higher dimension operators 
(D > 8) can be neglected safely. 

We adopt the naive dimensional regularization with the anticommuting 7 5 scheme, where 
there is no distinction between the first 4 dimensions and the remaining D — 4 dimensions. 
Since the bare effective Lagrangian contains the ultraviolet divergence which is induced 
by divergent subdiagrams, we give the renormalized results in the on-mass-shell scheme 



18| . Additional, we adopt the nonlinear gauge with £ = 1 for simplification [l9j. This 



special gauge-fixing term guarantees explicit electromagnetic gauge invariance throughout 
the calculation, not just at the end because the choice of gauge-fixing term eliminates the 
, yW ± G T vertex in the Lagrangian. 

Since the lepton EDM is an interesting topic in both theoretical and experimental aspects 
the current experimental upper limit on the electron EDM is 1.7 x 10~ 27 e • cm at 95% 
CL|21], and a future experiment with precision of 10 _29 e • cm is also proposed [221]. we as 
well present the lepton EDM by keeping all possible CP violating phases. Certainly, some 
diagrams in FigJI] have been discussed in Ref.jll] where the authors apply the projecting 
operator to get the lepton MDMs (Eq.8~Eq.10 in Ref.jjj]]). Nevertheless, the substantive 
corrections from several diagrams are ignored unreasonably (Fig. 5 in Ref. [III]). Additional, 
our formulae are new in their analytical forms. For the analysis on the electron EDM, we 
also include the contributions from the self energy diagrams (c) and (f) that are neglected 
in Ref . [3, Q] • Our result is more universal than that of Ref. {la. \h\ because we give up the 
assumption |mi j2 |, \/j, H \ ^> m w in concrete analysis. 

This paper is composed by the sections as follows. In section [Til we introduce the effective 
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Lagrangian method and our notations. Then we will demonstrate how to obtain the super- 
symmetric two-loop corrections to the lepton MDMs and EDMs. Section II III is devoted to 
the numerical analysis and discussion. In section IIVt we give our conclusion. Some tedious 
formulae are collected in appendix. 



II. NOTATIONS AND TWO-LOOP SUPERSYMMETRIC CORRECTIONS 

The lepton MDMs and EDMs can actually be expressed as the operators 

jC = 6 a la^l F 

MDM 4 m 1 » v ' 

^ edm = -\ d t l^lsl . (3) 

Here a^ u = z[/y M , 7„]/2, I denotes the lepton fermion, F^ v is the electromagnetic field strength, 
m l is the lepton mass and e represents the electric charge. Note that the lepton here is on- 
shell. 

In fact, it is convenient to get the corrections from loop diagrams to lepton MDMs and 
EDMs in terms of the effective Lagrangian method, if the masses of internal lines are much 
heavier than the external lepton mass. Assuming external leptons as well as photon are all 
off-shell, we expand the amplitude of the corresponding triangle diagrams according to the 
external momenta of leptons and photon. Then we can get all high dimension operators 
together with their coefficients. As discussed in the section HJ it is enough to retain only 
those dimension 6 operators in later calculations: 

of = ^r(<p)Vi, 

e 



°7 = ^ Wfj^ t i , 

^ eQ f m, - 
Of = ^lF.*^l, 



(4) 



with V = d +ieA and uj t = (1 =F7 5 )/2. When the equations of motion are applied to the 
incoming and outgoing leptons separately, only the operators Of 3 actually contribute to the 
MDMs and EDMs of leptons. We will only present the Wilson coefficients of the operators 
0^ 3 g in the effective Lagrangian in our following narration because of the reason mentioned 
above. We will adopt below a terminology where, for example, the "7/2-0" contribution means 
the sum of amplitude from those triangle diagrams (indeed three triangles bound together), 
in which a closed fermion (chargino/neutralino) loop is attached to the virtual Higgs and 
photon fields with a real photon attached in all possible ways to the internal lines. Because 
the sum of amplitude from those "triangle" diagrams corresponding to each "self-energy" 
obviously respects the Ward identity requested by QED gauge symmetry, we can calculate 
the contributions of all the " self-energies" separately. Taking the same steps which we did 

n n n 

in our earlier works [l_5|, lla, [17], we obtain the effective Lagrangian that originates from 
the self energy diagrams in FigHJ In the bare effective Lagrangian from the 'WW and 'ZZ' 
contributions, the ultraviolet divergence caused by divergent sub-diagrams can be subtracted 



safely in on- mass-shell scheme 



18] . Now, we present the effective Lagrangian corresponding 



to the diagrams in FigJT] respectively. 



A. The effective Lagrangian from 7/10 (7G0) sector 

As a closed chargino loop is attached to the virtual neutral Higgs and photon fields, a 
real photon can be emitted from either the virtual lepton or the virtual charginos in the self 
energy diagram. When a real photon is emitted from the virtual charginos, the corresponding 



"triangle" diagrams belong to the typical two- loop Bar-Zee-type diagrams 23]. Within the 
framework of minimal supersymmetric extension of the SM, the contributions from two- 
loop Bar-Zee-type diagrams to the EDMs of those light fermions are discussed extensively 
in literature 24(. When a real photon is attached to the internal standard fermion, the 
correction from corresponding triangle diagram to the effective Lagrangian is zero because 
of the Furry theorem, this point is also verified through a strict analysis. The corresponding 
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effective Lagrangian from this sector is written as 



with 



Clh 2 V2(4tt) 2 s2 A2 



x 

+iZ(H aa )(^-) 1/2 T 2 (x h ,x x ± )(o: - O;) 



(5) 



n a p = (Ul) a2 (U L ) lfj cos(3 + (Ul) al (U L ) 2p smP . 



(6) 



Where U L R denote the left- and right-mixing matrices of charginos, A denotes a energy scale 
to define Xi = mf/A 2 , respectively. The angle (3 is defined through the ratio between the 
vacuum expectation values of two Higgs doublets: tan (3 = v 2 jv\. We adopt the shortcut no- 
tations: c w = cos 9 w , s w = sin 6> w , where # w is the Weinberg angle. The concrete expressions 
of Ti )2 can be found in appendix. 

Accordingly, the lepton MDMs and EDMs from 7/io sector are written as 

,2 



7/1 



df 



(4tt)VA 2 1 ) l{Xh)X >cV xi> ' 

X 



v / 2(4tt) 4 s 2 A 2 

In the limit x ± x h , the above expressions can be simplified as 



(7) 



V2e 4 Q f 



m. 



■yh 



(4tt) 4 s 2 A 2 
v / 2(4vr) 4 s 2 A 2 



9fc(W««)(^) lim 



d 



x„, / x ±^ x ± dx 



v 2;,.. 7 



± \ 1/2 



lnx,. 



^ + lim — 

x . x ±-> x ± ox 



d 



-<P!{X ,,X . 



(8) 



Similarly, we can formulate the corrections from 7G0 sector to the effective Lagrangian 



as 



C "° 2 v / 2(47r) 2 s 2 A 2 ' 



X 

K(H aa )(^) 1/2 T 2 (x z ,x ± ,x J (O+ + O;) 

X„, 7 X a Xg \ 



^ QQ )(^) 1/2 r 1 (x z , x x - o;)\ . 



(9) 
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Correspondingly, the corrections to the lepton MDMs and EDMs from this sector are: 

,2 X 



7G 



x/2(4vr) 4 s 2 A 
In the limit x , ^> x , we have 



7G _ V2e"Q f mf , 
1 ~ (47r) 4 s 2 A 2 1 



.r 



± \ 1/2 



In a; , (9 
^ + lim — 

x , z i^.i ax 



(10) 



j7 g = e 5 Q/^ 



v/2(4tt) 4 s 2 A 2 



3(^0 



.7; 



± \ 1/2 



lim 



<9 



<9x 



-ipi[x 



(11) 



It should be emphasized that the corrections from this sector to the lepton EDMs are 
neglected in the analysis before Q, Q]. However, EqJTU] implies that the contributions from 
those diagrams to the lepton MDMs and EDMs can not be ignored generally 

Using the concrete expression of (fi(x, y) collected in appendix, one can verify easily that 
the corrections to the lepton MDMs and EDMs from the sectors are suppressed by the 
masses of charginos asm ± > m h , m z (a — 1, 2). 

B. The effective Lagrangian from Zho (ZGq) sector 



As a closed chargino loop is attached to the virtual Higgs and Z gauge boson fields, a 
real photon can be attached to either the virtual lepton or the virtual charginos in the self 
energy diagram. When a real photon is attached to the virtual lepton, the corresponding 
amplitude only modifies the Wilson coefficients of the operators Of in the effective La- 
grangian after the heavy freedoms are integrated out. In other words, this triangle diagram 
does not contribute to the lepton MDMs and EDMs. A real photon can be only attached to 
the virtual lepton as the closed loop is composed of neutralinos, the corresponding triangle 
diagram does not affect the theoretical predictions on the lepton MDMs and EDMs for the 
same reason. Considering the points above, we formulate the contributions from ZHq sector 
to the effective Lagrangian as 



x 



16\/2(47r)Vc 2 Q f A 2 



(T?~2Q f sl 



2(2 + lnx ± )g 0tl (x„x h ) 



+F 1 (x z ,x h ,x ± , x ± )Mn,J^ + HlJ%) (O: + Or) 



x 



+i(-^ L ) 1,2 \ -2(lnz ± - hxx ± )g 0tl (x z , x h ) + Fi{x a , x h , x ± ,x 

* Cl, / L X,-v X ,-) 7Cn ?i 



with 



^ = 25 a/3 cos2e w + (C/t) al ([/J li3 , 
e2 3 = 25 a/3 cos2e w + (C/t) ai ([/J i ^ 



(12) 



(13) 



where the concrete expressions of the functions ^ J (xi,x 2 ), F 1)2 (xi, x 2 , x 3 , x±) are listed in 
appendix. Additional, is the isospin of lepton, and Q f is the electric charge of lepton, 
respectively. Using EqJT2"l we get the corrections to the lepton MDMs and EDMs from Zh 
sector as 



Zh 



4 2 

em 



x 



4v / 2(4tt) 4 s 4 c 2 A 2 

V / WW 



(Tf-2Q f sl 



■J',. 



1/2 r 



2(2 + In x ± )g id {x x ,x h ) 



+F 1 (x z ,x h ,x x p(n j L a , + ^lCj > 



df h 



e 5 m, 



:{Tf-2Q f s 



1,1 4\V 2 



•X'.. 



2 (In a; ± - In a; ± )g Ax„,x h ] 



8v / 2(4tt) 4 s 4 c 2 A 2 

\ / w w 

+F 1 (x x , x hJ x , x ) + F 2 (x z , x k , x , x J] 3(7^ - Wj^j) • (14) 



The above equations contain the suppression factor 1— 4s 2 because Q f = —1 and - 
for charged leptons. In the limit x ± , x ± ^> x z , x h , Eq{14]can be approximated as 



-1/2 



Zh 



4 2 

em 



x 



(T?-2Q,sl){—f 2 [^^J 



4v / 2(4vr) 4 s 4 c 2 A 2 
2 2x ± g a (x ± , x 



d<pi 



OX , V ^ ' Xg ' 



e 5 m, 



8a/2(4tt)Vc 2 A 21 ' ^ J ^ x w J Kftc ± + dxJ KX & X *t } 



(15) 
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Similarly, the contributions from ZGq sector to the effective Lagrangian are 

„4 ( x 

-X 



-:(-^) 1/2 f-(2 + ln 



~^ x ± ) ~^ ^1 fax ' ~^x J x ± ' *^ x ± / 



lQV2(4n) 2 slclQ f A 2 { 
x%(h b i L +H) i R )(T z -2Q f s 2 )(0; -Of) 

+ (— ) [ (hix ± -In a; ± ) + Fi(x z ,x z ,x ± ,x J + F 2 (x z ,x z ,x ± ,x 

X v X z X a Xp Xa Xp Xp 



and the contributions to the lepton MDMs and EDMs are: 



a ZG = 



4 2 

em 



4v / 2(4tt)Vc 2 A 2 

V / w w 



{Tf-2Q f s 



x 

2\f X P 



x± x 1/2 r 2 



x„ 



(lnx ± — lnx 



x 



+Fi(x z ,x z ,x ± ,x J + F 2 (x z ,x z ,x ± ,x ± 



d? G = 



em 



-(T z -2Q f s 2 



X 



a;.. 



1/2 r 2 



8 v / 2(47t)Vc 2 A 2 

V 'WW 

+Fi(x.,x.,x^,x x± )]9(w <Ja ^+Wj a e. 



— (2 + lnx 



L x 



When x ± , x ± ^> x z , EqJTT] can be approached by 



^ZG 



4 2 

em 



x 



~4 v / 2(47r)Vc 2 A 2 

'WW 



(T/-2Q ;s! .)(-ll) 1/2 [( 



9^1 fy>l 



+ 



v <9x , <9x 



(x ± , X ± ) 

Xq. X 



d 



em 



8v / 2(4vr) 4 s 4 c 2 A 2 

\ 'WW 



(Tf-2g /S 2 )(^) 1/2 [|^(x ± ,xj 



• (1 + lnx )]Q(H H £ L ^ 



(16) 



(17) 



(18) 



C. The effective Lagrangian from 7Z sector 



When a closed chargino loop is attached to the virtual 7 and Z gauge bosons, the 
corresponding correction to the effective Lagrangian is very tedious. If we ignore the terms 
which are proportional to the suppression factor 1 — 4s 2 , the correction from this sector to 
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the effective Lagrangian is drastically simplified as 



' 7 Z 



8(4tt) 2 s 2 c 2 A 2 

\ / WW 



£ L -£ R ) lim T 3 (x ,x 



X 



Xa 



X 



(rf - QX)(° a ~ + O;) + Q f sl(0+ + 0+) 



+ 



(19) 



Correspondingly, the correction to the lepton MDMs from this sector is written as 



-yZ 

a = 



em 



■U L ~£ R ) lim T 3 (x ,x ± ,x . ) , 



4(47r) 4 S 2 C 2 A 2 ^' aa ' aa/ x ±^ x ± " " 

V / w W YT Yl 



(20) 



and the correction to the lepton EDMs is zero. In the limit x ± ^> x z , we can approximate 
the correction to the lepton MDMs from this sector as 



a] = 



4 2 

em 



4(4vr) 4 s 2 c 2 A 2 

V 'WW 



y ^ckck ^ act J 
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18x 



+ 



lnx ± — 2 lnx z 



3x 



Xa 



<)' 2 -: d(pi 
+ lim ( 2x + — ) [x ,,x 



Xa X 



= — = *i dx 2 . dx 



(21) 



D. The effective Lagrangian from WG^ sector 



As a closed chargino-neutralino loop is attached to the virtual W ± gauge boson and 
charged goldstone G T , the induced Lagrangian can be written as 



x 



CwG 16(4tt) Vc w Q,A 2 I ... 



4 \ V2 



f 3 (x w , Xw , x xl , x x± ) + QIC) o; 



+ i-f) 1/2 f,{x^x x% , x x± ) [(^c: - <?x) o fl - + ((<?% (c% 
-(e%(c*)LK] 

+ (^) 1/2 F 6 (x w ,a; w ,X x0a ,X x± )[(^Cf, - <?XK + ((S%(CT 

w P 

-(0%(c £ )J>. + ' 



(22) 
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with 

$t = ^{^(^(M^ + AT 2qCw ) - (U L ) 10 Af 3a c w } , 

Gf a =-cosf3{±(Ul) 02 (^ lSw +Ml 2 c v ) - (UD^AT^} . (23) 

Here, the 4x4 matrix M denotes the mixing matrix of the four neutralinos Xa i a = 
1, 2, 3, 4). 

The corresponding corrections to the lepton MDMs and EDMs are respectively expressed 

as 



w 

X x° \ V 2 



F 5 (x ,x ,x n ,x + )K(£ L C L - G R C R 



/3a a/3 {3 a a/3 



2 f 6 (x ,x ,x n , x + w(g L c R -g r c l )\ 

+ (^) 1/2 F 4 (x w ,x w ,x x J^X+^C 



Using the asymptotic formulations of form factors 245 6,7 collected in appendix, we can 
simplify the expressions of EqfSH in the limit x , x ± 3> x w . 

The contributions from those above sectors to effective Lagrangian do not contain ul- 
traviolet divergence. In the pieces discussed below, the coefficients of high dimensional 
operators in effective Lagrangian contain ultraviolet divergence that is caused by the di- 
vergent subdiagrams. In order to obtain physical predictions of lepton MDMs and EDMs, 

13 



it is necessary to adopt a concrete renormalization scheme removing the ultraviolet diver- 
gence. In literature, the on-shell renormalization scheme is adopted frequently to subtract 



the ultraviolet divergence which appears in the radiative electroweak corrections [181 ] . As 
an over-subtract scheme, the counter terms include some finite terms which originate from 
those renormalization conditions in the on-shell scheme beside the ultraviolet divergence to 
cancel the corresponding ultraviolet divergence in the bare Lagrangian. In the concrete cal- 
culation performed here, we apply this scheme to subtract the ultraviolet divergence caused 
by the divergent subdiagrams. 



E. The effective Lagrangian from the ZZ sector 

The self energy of Z gauge boson composed of a closed chargino loop induces the 
ultraviolet divergence in the Wilson coefficients of effective Lagrangian. Generally, the un- 
renormalized self energy of the weak gauge boson Z can be written as 

X? (p) = k 2 Alg, u + (A\ + ^ A*) V - P»P») + {Bl + £BZ)p&, . (25) 
Correspondingly, the counter terms are given as 

^(p) = ~M + m l SZ *)9^ ~ - P&v) • (26) 

The renormalized self energy is given by 

£>) = £>) + 5t c (p). (27) 
For on-shell external gauge boson Z, we have [11 



o 



lim -jJ— jl« (p)^(p) = c (p) , (28) 

z ± z 

where e(p) is the polarization vector of Z gauge boson. From Eq. (128]) . we get the counter 
terms 

2 

Til 

5Z z =A z 1 + 1 ^A z 2 = Al + xA z 2 , 

dm] = A Z Q A 2 - m 2 5Z z . (29) 
14 



-*s zc 

I '"' I 




FIG. 2: The counter term diagram to cancel the ultraviolet caused by the self energy of Z boson. 



Accordingly, the effective Lagrangian originating from the counter term diagram (Figf2]) 
can be formulated as 



5C 



c 



2£ r^i ±£ ) 

12(47r)Vc 4 A 2Vi "~ Rj (l-e) 2 



r (47TX p 



e x z 



5[x ± -\- x ± ) Q 2 ,i( x ±-> x ±) 



12a; 2 



l£ I Xa X f3 _|_ 

X 2 



iC _j_ I tjC _|_ 

X a X^ 

H ^ — - mx t 



+2(x , x 



Xa X ^ 



12a; ar 

z z 

1 #1,1 ± 5 X ±) 

L Xa X^ 

~ 2 



+ 



In a^ 



ex 

2 A f/n+ 



X, 



12a; 2 



T 2 (l + e ' 



r (47rx R 



,2e- 



4(4vr) 2 s 4 c 4 A 2V Ry (l-e) 2 

\ / WW v ' 

Q2 1 -j- ; ^ -j- ) X _j_ ~T~ X _j_ y 

**-( — h In a;, — lnx ) H 

v 2 1 ° J Ax, 



1 



£ X* 

Z 

1 I 3j 



In a; T 



a; z 



a 



x- 



+2(x x x± ) 1/2 (CC+C< 



^1,1 ± ) ^ ± ) 



ear 



ar 



4c^-^ + !!?kw-«)w^-- (30, 



a;" 



ar 



Here, £ = 2 — D/2 with .D representing the time-space dimension, and x a = A /A (A RE 
denotes the renormalization scale). 

As a result of the preparation mentioned above, we can add the contributions from the 
counter term diagram to cancel the corresponding ultraviolet divergence in bare effective 
Lagrangian. The resulted theoretical predictions on the lepton MDMs and EDMs are re- 
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spectively written as 



n ZZ — 



4 2 

era 

i 



(4tt) 4 s 4 c 4 A 2 

\ / WW 



(iCi 2 + iCi 2 )[( T /-e,<) 2 +^ 



-± (r 5 (x, , x ± , x ± ) + "° 2 - Ik,)+jTi(i„i„iJ 

+1OCI 2 - ICI 2 ) [ W - <?X) 2 - T ^\i • \j ) 



X 



4Q, x z 7Q f 



T 7 (x z ,x x ) + -Zi.ln^- 

"^CK -*/3 OX X-r. OX 



-(lez+ie/KC^-e,- 

z 

3 



g 2 



2a; 2 



2x 2 



{q 21 {x ±i% ±) ^ ±# ±f?o,i( 3 ' ± ' ±)) 



x 



(4tt) 4 s 
d 2 



4 

w w 



<9 2 



dxj)x ± dx z dx ± 
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Because a real photon can not be attached to the internal closed neutralino loop, the 
corresponding effective Lagrangian only contains the corrections to the lepton MDMs: 
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X 

"2(^ o "I" o ) (2 — l n ) ~~ f?2,i q j x v o ) Qo,l \ x „o j ) 



-«?,*(«<x( r ; - o^x**^) 1 " 2 'y' 1 ' 1 1 (32) 



with 



Va(3 = A/J^ , 

^=^3^,(^/3 = 1, ■••,4). (33) 

We can also simplify Eq. ([3"Tj) and Eq. (j52|) using the asymptotic expressions of T 8 ~ T 13 in 
the limit x ± , x ± , x , x ^ x z . The concrete expressions of T 8 ~ T 13 can be found in 
appendix. 



X a Xp X a 



F. The effective Lagrangian from the WW sector 

Similarly, the self energy of W gauge boson composed of a closed chargino-neutralino 
loop induces the ultraviolet divergence in the Wilson coefficients of effective Lagrangian. 
Accordingly, the unrenormalized W self energy is expressed as 

£W(p) = A 2 A-g, u + (AY + ^A*) {p 2 g, u - p,p v ) + (flf + ^B-)p, Pl/ . (34) 

The corresponding counter terms are given as 

*V°(P) = ~( 6m l + < 5Z ^ - SZrfgp, - p, Vv ) . (35) 
The renormalized self energy is given by 

S» = S» + E^(p) (36) 
For on-shell external gauge boson W ± , we have [3] 

Sj(p)e"(p) 2 2 = 0, 

1 w 

2 lim 2 ^T^(p) e » = 6,(p), (37) 

W 1 w 
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FIG. 3: The counter term diagram to cancel the ultraviolet caused by the self energy of W boson 
and electroweak radiative corrections to ■jW + W~ vertex. 

where e{p) is the polarization vector of W gauge boson. Inserting Eq. and Eq. (|35|) 
into Eq. (1371) . we derive the counter terms for the W self energy as 

2 

Try 

6ml = AqA 2 - m 2 JZ w . (38) 

Differing from the analysis in the ZZ sector, we should derive the counter term for the vertex 
r )W + W~ here since the corresponding coupling is not zero at tree level. In the nonlinear 
gauge with £ = 1, the counter term for the vertex r jW + W~ is 

iSC lW + w - = ie ■ SZ w [g^(h - k 2 ) p + g up (k 2 - h)^ + g Pfl (k 3 - fc^] , (39) 

where ki (i — 1, 2, 3) denote the injection momenta of W ± and photon, and /i, u, p denote 
the corresponding Lorentz indices respectively. 
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We present the counter term diagrams to cancel the ultraviolet divergence contained in 
the bare effective Lagrangian from WW sector in FigfSl and we can verify that the sum of 
corresponding amplitude satisfies the Ward identity required by the QED gauge invariance 
obviously. Accordingly, the effective Lagrangian originating from the counter term diagrams 
can be written as 
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Finally, we get the renormalized effective Lagrangian from the WW sector: 
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(41) 



Correspondingly, the resulted lepton MDMs and EDMs are respectively formulated as 
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III. NUMERICAL RESULTS AND DISCUSSION 

With the theoretical formulae derived in previous section, we numerically analyze the 
dependence of the muon MDM and the electron EDM on the supersymmetric parameters in 
the split scenario here. In particular, we will present the dependence of the muon MDM and 
the electron EDM on the supersymmetric CP phases in some detail. Within three standard 
error deviations, the present experimental data can tolerate new physics corrections to the 
muon MDM as -10 x 10~ 10 < Aa M < 52 x 10~ 10 . Since the neutralinos xl (« = 1, 2, 3, 4) 
and charginos xt i a = 1j 2) appear as the internal intermediate particles in the two-loop 
diagrams which are investigated in this work, the corrections of these diagrams will be 
suppressed strongly when the masses of neutralinos and charginos are much higher than 
the electroweak scale Uj . To investigate if those diagrams can result in concrete corrections 
to the muon MDM and electron EDM, we choose a suitable supersymmetric parameter 
region where the masses of neutralinos and charginos are lying in the range M x < 500 GeV. 
Without losing too much generality, we assume the supersymmetric parameters satisfying 
| m-i | = |m 2 | in this work. In split SUSY, only the CP violating phases arg(/i ff ), arg(mi), 
and arg(m2) have substantive contributions to the lepton MDMs and EDMs. As for other 
CP violating phases, either they do not contribute to the lepton MDMs and EDMs or their 
corrections can be neglected safely. Moreover, the existence of a CP-even SM like Higgs with 
mass around 100 — 250 GeV sets a strong constraint on the parameter space of the employed 
b address this problem, they argue that a fine tuning in the Higgs potential is 



model. 

required [25}]. Admitting this fine tuning among high energy scale parameters, one no longer 
worries about the constraint from Higgs sector. On the other hand, the CP violation would 
cause changes to the neutral-Higgs-quark coupling, the neutral Higgs-gauge-boson coupling 
and the self-coupling of Higgs boson. The present experimental lower bound on the mass 
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FIG. 4: The supersymmetric corrections to the electron EDM d e and muon MDM vary with 
the CP violating phase 4> CP when \fj, H \ = \nii\ = \m2\ = 200 GeV and tan/3 = 5. Where the dash 



lines stand for the corrections with 



arg(mi) and arg(m,2) = avg(fi H ) = 0, the dot lines stand 



for the corrections with 



the corrections with 



CP 



arg(m2) and arg(mi) = arg(/x H ) = 0, and the solid lines stand for 



arg(n H ) and arg(mi) = arg(m2) = 0, respectively. 



of the lightest Higgs bosons is relaxed to 60 GeV 26j. The scope of Higgs mass is chosen 
around 60 — 250 GeV in our numerical analysis because of the above reasons. In fact, we 
find that the lepton MDMs and EDMs weakly depend on the mass of lightest Higgs by 
scanning the parameter space. In the following discussion, we choose the mass of light Higgs 
as m h = 120 GeV. 

Taking = \mi\ = \m2\ = 200 GeV and tan/3 = 5, we plot the electron EDM d e and 
muon MDM a M versus the CP phases CP = arg(mi), arg(m2), arg( y u J/ ) separately in FigJH 
When arg(m 2 ) = arg(/i H ) = 0, there is cancellation among the dominant contributions to 

that originate from the "7Z" and "ZZ" sectors respectively. As 4> CP = aig(fi H ) and 
arg(mx) = arg(m 2 ) = 0, the absolute values of supersymmetric corrections to the electron 
EDM d e (solid line in FigJUa)) exceed 2.5 x 10~ 28 e-cm at the largest CP violation arg(/i^ 
7r/2, 37r/2, which is well below the present experimental upper limit 1.7 x 10~ 27 e • cm 



2l|, 
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but large enough to be detected in next generation experiments. Correspondingly, the muon 
MDM a M depends on the CP violating phase CP = arg(/x H ) (solid line in FigJU^b)) strongly, 
the super symmetric corrections to exceed 3 x 10~ n at the largest CP violation (arg(/i ff ) = 
7r/2, 3tc/2) since the cancellation mentioned above dissolves now. For the same reason, the 
super symmetric correction to surpasses 7 x 10~ n at the CP conservation of arg(yU H ) = ir. 
As cf) CP = arg(m 2 ) and arg(mi) = arg(^ H ) = 0, the absolute values of supersymmetric 
corrections to the electron EDM d e (solid line in FigJD^a)) exceed 2 x 10 -28 e-cm at the largest 
CP violation arg(/! H ) = n/2, 3n/2, which is expected to be observed in next generation 



experiments with the sensitivity 10~ 29 e • cm 22|. Correspondingly, the muon MDM 
depends on the CP violating phase (p CP = arg(m 2 ) (dot line in Figfjjb)) sensitively, the 
supersymmetric corrections to a M are about 3 x 10 -11 at the largest CP violation (arg(/i ff ) = 
7r/2, 37r/2) because the cancellation existing as arg(m 2 ) = arg(/i H ) = vanishes here. When 
<ft CP = arg(mi) and arg(m 2 ) = arg(/z H ) = 0, the supersymmetric correction to the electron 
EDM d e (dash line in FigJU^a)) is below 1 x 10 -30 e • cm, and very difficult to be detected 
in near future. Corresponding to small theoretical prediction on the electron EDM d e , the 
muon MDM varies with the CP phase CP = arg(mx) (dash line FigJU^b)) very mildly. 

Taking \fi H \ = \mi\ = |m 2 | = 200 GeV and tan/3 = 50, we plot the electron EDM d e and 
muon MDM a M versus the CP phases CP = arg(mi), arg(m 2 ), arg(/i H ) separately in FigJHl 
As CP = arg(/i H ) and arg(mi) = arg(m 2 ) = 0, the absolute values of supersymmetric 
corrections to the electron EDM d e (solid line in Figjf^a)) reach 3 x 10 -28 e • cm at the 
largest CP violation arg(yU ff ) = tt/2, 3tt/2, which exceeds the precision of next generation 
experiments[22]. Correspondingly, the muon MDM depends on the CP violating phase 
4> CP = arg(/i ff ) (solid line in Figj5]^b)) strongly, the supersymmetric corrections to a M are 
about 3 x 10 -11 at the largest CP violation (arg(/i H ) = tt/2, 37r/2) since the cancellation 
appearing at arg(m 2 ) = arg(/i ff ) = dissolves here. As (fi CP = arg(m 2 ) and arg(mi) = 
arg(yU H ) = 0, the absolute values of supersymmetric corrections to the electron EDM d e (solid 
line in Figj5](a)) exceed 2.5 x 10~ 28 e ■ cm at the largest CP violation arg(/i H ) = n/2, 37r/2, 
which is expected to be observed in next generation experiments with the sensitivity 10 -29 e- 
cm. Correspondingly, the muon MDM a M depends on the CP violating phase CP = arg(m 2 ) 
(dot line in Figj5]Jb)) steeply, the supersymmetric corrections to a M are about 4 x 10~ n at the 
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FIG. 5: The supersymmetric corrections to the electron EDM d e and muon MDM vary with 
the CP violating phase <p CP when \/j, H \ = \m\\ = \m2\ = 200 GeV and tan/3 = 50. Where the dash 
lines stand for the corrections with <f> CP = arg(mi) and arg(m2) = arg(/x H ) = 0, the dot lines stand 
for the corrections with (p CP = arg(rri2) and arg(mi) = arg(/x H ) = 0, and the solid lines stand for 
the corrections with <p CP = arg(/z H ) and arg(mi) = arg(m2) = 0, separately. 

largest CP violation (arg(/x H ) = tt/2, 37r/2) because the cancellation existing at arg(m2) = 
axg(fi H ) = disappears presently. When CP = arg(mi) and arg(m 2 ) = arg(/i H ) = 0, 
the supersymmetric correction to the electron EDM d e (dash line in FigJ3](a)) is less than 
1 x 10~ 30 e • cm, and very difficult to be detected in near future. Corresponding to small 
theoretical prediction on the electron EDM d e , the muon MDM varies with the CP phase 
CP = arg(mi) (dash line FigJS^b)) very slowly. 

Note that the theoretical predictions on the electron EDM d e and muon MDM a M are not 
enhanced by large tan j3 here, this point can be understood as follows. The tan f3 enhanced 
couplings are only contained in the interactions among the heavy Higgs and down-type 
fermions (sfermions). Those heavy Higgs fields include the neutral CP-odd Higgs, the neutral 
heavy CP-even Higgs, as well as the charged Higgs. However, those particles are decoupled 
from the low energy theory because they are super-heavy under the split assumption. In 
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FIG. 6: The supersymmetric corrections to the electron EDM d e and muon MDM vary with 
the CP violating phase 4> CP when \fj, H \ = \nii\ = |m 2 | = 300 GeV and tan/3 = 5. Where the dash 
lines stand for the corrections with </> CP = arg(mi) and arg(m 2 ) = arg(/x ff ) = 0, the dot lines stand 
for the corrections with (p CP = arg(m2) and arg(mi) = arg(/x H ) = 0, and the solid lines stand for 
the corrections with CP = arg(^ H ) and arg(mi) = arg(m2) = 0, separately. 

other words, our theoretical predictions are not enhanced by large tan /3 since we ignore the 
corrections from those heavy Higgs fields. 

Taking \fi H \ = \mx\ = |m 2 | = 300 GeV and tan/3 = 5, we plot the electron EDM d e 
and muon MDM a M versus the CP phases CP = arg(mi), arg(m 2 ), arg(/i H ) separately in 
FigEl As arg(m 2 ) = arg(/i £r ) = 0, a cancellation exists among the dominant supersymmetric 
contributions to a M that originate from the "7Z" and " ZZ n sectors respectively. When </> CP = 
arg(yU £r ) and arg(mi) = arg(m 2 ) = 0, the absolute values of supersymmetric corrections to 
the electron EDM d e (solid line in Figj6]^a)) approach 2 x 10~ 28 e • cm at the largest CP 
violation arg(/! H ) = it/ 2, 3ir/2, which is well below the present experimental upper limit 
1.7 x 10~ 27 e • cm, but large enough to be detected in next generation experiments with 
the precision of 10 x 10~ 29 e • cm. Correspondingly, the muon MDM depends on the 
CP violating phase CP = arg(/i H ) (solid line in Figj6](b)) strongly, the supersymmetric 
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corrections to a M exceed 2 x 1CT 11 at the largest CP violation (arg(/i H ) = ir/2, 3n/2) since 
the cancellation mentioned above dissolves here. For the same reason, the supersymmetric 
correction to surpasses 5.5 x 10~ n at the CP conservation of arg(/! ff ) = n. As <f) CP = 
arg(m 2 ) and arg(mi) = arg(/i H ) = 0, the absolute values of supersymmetric corrections to 
the electron EDM d e (solid line in Figj6]Ja)) exceed 2 x 10~ 28 e-cm at the largest CP violation 
arg(yU £r ) = 7r/2, 3tt/2, which is expected to be observed in next generation experiments. 
Correspondingly, the muon MDM a M depends on the CP violating phase CP = arg(m 2 ) 
(dot line in FigJI^b)) steeply, the supersymmetric corrections to a M are about 2 x 10~ n 
at the largest CP violation (arg(/x H ) = tc/2, 3tc/2) because the cancellation appearing at 
arg(m 2 ) = Big(fi H ) = vanishes now. When CP = arg(mi) and arg(m 2 ) = arg(/i H ) = 0, 
the supersymmetric correction to the electron EDM d e (dash line in FigEfa)) is below 

1 x 10~ 30 e • cm, and very difficult to be detected in near future. Corresponding to small 
theoretical prediction on the electron EDM d e , the muon MDM a M varies with the CP phase 
CP = arg(mi) (dash line FigJSfb)) very slowly. 

Taking \fi H \ = \rri\\ = |m 2 | = 300 GeV and tan/5 = 50, we plot the electron EDM d e 
and muon MDM a M versus the CP phases CP = arg(mi), arg(m 2 ), arg(yU H ) separately in 
Figl3 As CP = arg(/x H ) and arg(mi) = arg(m 2 ) = 0, the absolute values of supersymmetric 
corrections to the electron EDM d e (solid line in FigH^a)) reach 2.3 x 10 -28 e ■ cm at the 
largest CP violation arg(/i H ) = tt/2, 3tt/2, which exceeds the precision of next generation 
experiments. Correspondingly, the muon MDM a M depends on the CP violating phase <j) CP = 
arg(/i H ) (solid line in Fig 13(b)) strongly, the supersymmetric corrections to a M are about 

2 x 10~ n at the largest CP violation (arg(/i H ) — tt/2, 3tt/2) since the cancellation existing 
at arg(m 2 ) = arg(yU H ) = dissolves here. As CP = arg(m 2 ) and arg(mi) = arg(yU ff ) = 0, the 
absolute values of supersymmetric corrections to the electron EDM d e (solid line in FigJT^a)) 
exceed 2.2 x 10~ 28 e-cm at the largest CP violation arg(yU H ) = n/2, 3tc/2, which is expected 
to be observed in next generation experiments. Correspondingly, the muon MDM a M depends 
on the CP violating phase <f) CP = arg(m 2 ) (dot line in FigJT^b)) steeply, the supersymmetric 
corrections to are about 2 x 10~ n at the largest CP violation (arg(// H ) = tt/2, 3tt/2) 
because the cancellation mentioned above vanishes also now. When CP = arg(mi) and 
arg(m 2 ) = aig(fi H ) = 0, the supersymmetric correction to the electron EDM d e (dash line 
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FIG. 7: The supersymmetric corrections to the electron EDM d e and muon MDM vary with 
the CP violating phase <p CP when \/j, H \ = \m\\ = \m2\ = 300 GeV and tan/3 = 50. Where the dash 
lines stand for the corrections with <f> CP = arg(mi) and arg(m2) = arg(/x H ) = 0, the dot lines stand 
for the corrections with (p CP = arg(rri2) and arg(mi) = arg(/x H ) = 0, and the solid lines stand for 
the corrections with <f) CP = arg(/z H ) and arg(mi) = arg(m2) = 0, respectively. 

in FigJT^a)) is less than 1 x 10~ 31 e • cm, and very difficult to be detected in near future. 
Corresponding to small theoretical prediction on the electron EDM d e , the muon MDM a M 
varies with the CP phase 4> CP = arg(mi) (dash line FigJ7[b)) very mildly. 

In the numerical analysis presented above, the assumption \fi H \ = \mi\ = \m2\ is taken 
for simplicity. This assumption on parameter space induces very specific mixing patterns 
in chargino and neutralino sectors respectively. In order to investigate the supersymmetric 
corrections to and d e without the assumption, we plot a M and d e varying with the energy 
scale of new physics A NP when arg(/i H ) = arg(mi) = 0, arg(m 2 ) = vr/2 in FigJHl Because the 
supersymmetric corrections to and d e depend on tan (3 mildly, we choose a middle value of 
tan p = 20. In FigJS](a) and FigJS](b), the solid lines stand for the supersymmetric corrections 
with \/j, H \ — | mi | = 1 777.2 1 = A NP , the dot lines stand for the supersymmetric corrections with 
\fi H \ = 200 GeV, | mi | = \m2\ = A NP , and the dash lines stand for the supersymmetric 
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FIG. 8: The supersymmetric corrections to the electron EDM d e and muon MDM vary with the 
energy scale A when tan/3 = 20, and aig(fi H ) = arg(mi) = 0, arg(m2) = tt/2. Where the dot lines 
stand for the corrections with \fi H \ = 200 GeV, \m\\ = \m2\ = A NP , the dash lines stand for the 
corrections with \mi\ = \m2\ = 200 GeV, \n H \ = A NP , and the solid lines stand for the corrections 
with \fj, H \ = | mi | = \rri2\ = A NP , respectively. 

corrections with \m\\ = |m 2 | = 200 GeV, \n H \ = A NP , respectively. In FigJH^a), the 
resonance around A NP = 200 GeV is arisen by the intervention between the standard and 
supersymmetric fields. As A NP < 350 GeV, the difference between the theoretical predictions 
on d e with and without \[i H \ = \m>i \ = \m2\ is not very obvious. With the increasing of A NP , 
the suppression of supersymmetric correction to d e without the assumption is more stronger 
than that with the assumption. A similar case exists in the 2-loop electroweak correction to 
a^, the suppression of supersymmetric correction to a M without the assumption is stronger 
than that with the assumption when the energy scale A NP increases. 
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IV. CONCLUSIONS 



In this work, we analyzed the two-loop supersymmetric corrections to the muon MDM 
and electron EDM by the effective Lagrangian method in split scenarios. In the concrete 
calculation, we keep all dimension 6 operators. The ultraviolet divergence caused by the 
divergent sub-diagrams is removed in the on-shell renormalization schemes. After applying 
the equations of motion to the external leptons, we derive the muon MDM and the electron 
EDM. Numerically, we analyze the dependence of the muon MDM a M as well as the electron 
EDM d e on supersymmetric CP violating phases. Adopting our assumptions on parameter 
space of the split supersymmetry, we find that the correction from those two-loop diagrams 
to a M is below 10~ 10 roughly for CP conservation, which is less than the present experimental 
precision in magnitude. In other words, the present experimental data do not put a very 
restrictive bound on parameter space of split supersymmetry. Additional, the contribution 
to dfj, from this sector is sizable enough to be experimentally detected with the experimental 
precision of near future. 

Acknowledgments 

The work has been supported by the National Natural Science Foundation of China 
(NNSFC) with No. 10675027. One of us (TFF) also acknowledges the support from the 
ABRL Grant No. R14-2003-012-01001-0 of Korea at the early stage of this work. 



APPENDIX A: THE FUNCTIONS 

We list the tedious expressions of the functions adopted in the text 
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2 2 



+( 



a 3 



X1X3 x\ 

+ n o„ + 



(lnx 2 — lnx 3 ) 



dx\dx 3 dxidx 3 dx\dx 2 dx\dx 2 dx 3 



$(xi,x 2 ,x 3 ) 



3^2 - ^3 
Xi 

9 3 $ 



($(xi, x 2 , x 3 ) - ip (x 2 , x 3 )) 



r 8(ll ,, 2 , l3 ) = -4( s ^- + 



<9 3 $ 



dx\dx 3 dx\dx 2 
/ <9 3 <9 3 
V 9xi<9x| dx\dx 2 

-®(X!,X 2 ,X 3 ) 



X!,X 2 ,X 3 ) + 



+ -^(2 + lnx 2 ) 



X1X3 xf 

' 3 ($(xi, x 2 , x 3 ) - tp Q (x 2 , x 3 )) 



Xi 



r9(I1 , I2 , I3) .l ln , s _^ ( » ( ^ 2 ,, 3) _ s(I2 ,, 3)) 



X 2 

<9 2 /, $(xi,x 2 ,x 3 ) - <po(x 2 ,x 3 ) ^ f ,\ 

x 2 - x 3 ) $(Xi,X 2 ,X 3 ) J 



dxidx 



+-( 



X! V^X 3 ^Xi 

i 17x 2 

xf 

10(x 2 + x 3 ) 



(xi,x 2 ,x 3 ) + 



Xi 

4x3 <9 2 $ 



(xi,x 2 ,x 3 ) , 



. 26 17x 2 29x 3 10 . 
T w {x 1 ,x 2 ,x 3 ) = 1 2— H 2— + -^Q 21 {x 2 ,x 3 ) 

X ^ X ^ X ^ X ^ 



Xi dxidx 3 

16(x 2 - x 3 ) 2 



x 3 



X 2 



+ 
+ 
+ 



4 + 



61nx 3 
mxi h 

Xi 

16 (x 2 — x 3 )n X3 lnx 3 



14- 



16(x 2 - £3)1 x 2 lnx 2 



Xi 



x 2 



Xi 



(x 2 - x 3 ) 2 - x\ 



dxf 



Xi,X 2 ,X 3 ) 



5xi + 6x 2 + 



3(x 2 - x 3 ) ; 

Xi 



9 3 $ 



<9xf 



T (xi,x 2 ,x 3 ) 



9(x 2 -x 3 ) 2 t 6x 2 { 3x 3 -i<9 2 $ 

dx\ 



2 
xf 



+ — + 

Xi Xi 



(xi,x 2 ,x 3 ) 
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+ 
+ 



I2x 2 6x 3 18(x 2 -x 3 ) 2 i <9$ 



x\ 

12x2 6x 3 



x\ 



+ 



-(xi,X 2 ,X 3 ) 



18(x 2 -x 3 ) 2 i 



xf 



xf 



x\ 



dx 1 

($(x 1 ,x 2 ,x 3 ) - ipo(x 2 ,x 3 



2x\{x2 - X 3 ) r<9 3 $ 



xf 



dx\ 



(xi,X 2 ,X 3 ) 



dx 3 



2 ,x 3 ) 



+ 



x\ 



94-1 r 9 2 $ 

<9x 2 



XT 



Xi,X 2 ,X 3 ) 



(9X3 



(x 2 ,x 3 ) 



3x n 



xf 



+ 



9X/3 18x 3 (x 2 - x 3 ) 



x? 



+ 



xf 



<9$ 



<9x 3 



(xi,x 2 ,x 3 ) 



dx„ 

-6X3 (x 2 + X 3 - Xi) 



<9xf<9x 3 



-2x1 ( 



1 + 



x 2 - x 3 - 

Xi ' 



9 4 $ 

dx\dx\ 

<9 4 $ 



(Xi,X 2 ,X 3 ) 

: (xi,x 2 ,x 3 ) 



+ 
+ 



3xi — 3x 2 — I8X3 



dx\dx\ 

9x 3 (x 2 -x 3 )-i <9 3 $ 



- 21x 3 h 

Xi X\ 



Xi 

9x?i <9 3 $ 



<9xi<9x 2 



dx\dxz 
(xi,x 2 ,x 3 ) 



(xi,x 2 ,x 3 ) 



6 



12x 2 6x3 18x 3 (x 2 - x 3 ) 



T n (xi,x 2 ,x 3 ) 



<9 2 $ 

dxidx 3 

21nx 3 4(x 2 -x 3 ) 4(x 2 lnx 2 - x 3 lnx 3 ) 



+ 

Xi Xi 



(xi,x 2 ,x 3 ) , 



Xi 



xf 



4(x 2 - x 3 ) / . \ 4(x 2 - x 3 ) <9$ 

-3 — ($(xi,x 2 ,x 3 )-<Mx 2 ,x 3 )J + — - 2 — t^t(xi, 

2(x 2 - x 3 )n d 2 & 



x 2 



■(1+ 



\ u v 2x 3 / 9$ 

J ^(xi,x 2 ,x 3 ) - —(—( Xl ,x 2 ,x 3 ) 



Xi 



<9x 3 
<9 2 $ 



(x 2 ,x 3 )) + 



<9xf 

x 3 (x 2 - x 3 ) /<9 2 $ 



xf ^<9x 2 



xf 



^2 ( X 1,X 2 ,X 3 ) 



dx 2 



(x 2 ,x 3 )) 



<9xi<9x 3 

+ (x 2 + X 3 - Xi 



(xi,x 2 ,x 3 ) - X 3 (l + 

9 3 $ 



x 2 - x 3 



dx\dx 3 



Xi 

(xi,x 2 ,x 3 ) , 



<9 3 $ 

^Xl^Xg 



(xi,x 2 ,x 3 ) 
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Ti 2 (xi,x 2 ,x 3 ) = g 

Xf X1X3 



20, 181nx 3 20 

-2 lnxi 5 2-^ii(x 2 ,x 3 ) 



3 (<P(Xi,X 2 ,X 3 ) - ipo{X2,X 3 ) J + ^2 — (Xi,X 2 ,X 3 ) 

x 1 



xf dx 1 
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6 <9 2 $ 



<9 3 $ 



<9 4 $ 



<9x 3 ' ' <9x 3 
3(X2-2X3)/9 2 <I> ,\ 

15 to (^1,^2, ar 3 ) - -^(x^J 



x 3 (x 2 - x 3 ) /<9 3 $ <9 3 <£o 



x 2 



(, x 2 -x 3 ^ <9 4 $ 
"M 1 ; J nT^Is^l'^,^) 



xi ' <9xi<9x 3 
2(x 2 -x 3 )x <9 2 $ 



_6_/ -X' 3 K 
Xi ^ X\ 



dx±dx 3 
x 2 -2x 3 \ <9 3 $ 



3(l- X2 ~ ZX3 ) 



(Xl,X 2 ,X 3 ) 

x 2 — x 3 \ <9 3 $ 



6(2 - 



X! ' dx\dx 3 ^ 



x\ ) dx\dx\ 

<9 4 $ <9 4 $ 
+3(x 2 -x 3 - xi) (xi,x 2 ,x 3 ) -6 (xi,x 2 ,x 3 ) , 



(xi,x 2 ,x 3 ) 



^13(^1, x 2 , x 3 ) 



xix 3 x 2 ^<9x 3 



<9x 3 <9x 3 

dx 3 



<9x 2 <9x 2 



2 <9 2 $ 



(xi,x 2 ,x 3 ) - — -(x 2 ,x 3 )) - 

<9x 3 7 xi <9xi<9x 3 

x 2 -x 3 /<9 2 $ <9Vo, x\ 

— (-^2 (^1^2, ^3) - ^2-(x 2 ,X 3 )J 



(Xi,X 2 ,X 3 ) 



-0 



X! - L/x 3 

x 2 — x 3 \ <9 3< £> 



Fi(xi,x 2 ,x 3 ,x 4 ) 



xi 7 <9xi<9x 2 
J d_ 

XiX 2 9X4 

J j- 8 r/ 



(xi,x 2 ,x 3 ) - 2- 



<9 3 $ 



Xi — X 2 ^ (9X4 

9 



<9x 2 <9x 3 

x 3 - x 4 )v9o)(x 3 ,x 4 ) 

x 3 - X 4 \ -I 

J<Pj (Xi,X 3 ,X 4 ) 



(xi,x 2 ,x 3 ) , 



Xi 



(l + ^^ i )<f](x 2 ,x 3 ,x 4 )}, 



^2 



F 2 (xi,x 2 ,x 3 ,x 4 ) 



(9x 4 

J d_ 

xix 2 (9x 4 

J ( d r/ _ 

Xi — X 2 (9X4 ^ 



(x 3 - x 4 )y9 )(x 3 ,x 4 ) 

X 3 -X 4 x -I 

J<£>J (Xi,X 3 ,X 4 ) 



Xi 



9 



(9x4 



(l - — — — ) ®\(x 2 , x 3 , x 4 )} , 
^ x 2 ' J J 



F 3 (xi,x 2 ,x 3 ,x 4 ) = 2(lnx 4 - l)^ 01 (xi,x 2 ) - 

XiX 2 

xix 2 + 2(xi + x 2 )(x 3 - x 4 ) 



6(x 3 — x 4 ) 6(x 3 lnx 3 — X4I11X4) 



XiX 2 



+ 



X]^X 2 



x 3 - 3x 4 d(f 
<Po{xa,X4) ^—(^3,^4) 
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x 4 (x 3 - x 4 ) gVo 



X\X 2 



(x 3 ,x 4 ) 



_d_ 

dxi 



+ X 4 



J!! 

<9xj 



f2 (xi, X2; x 3 , x 4 ) 



M ^4' 

9 d 2 
+ (l - (^3 -3x 4 )- x 4 (x 3 - x 4 )—^jVL_ 1 (xi,x 2 ;x 3 ,x 4 ) 



-( 



<9 



+ 



\2r 



^(rci,^; a?3, ^4) + (£3 - £4)^0 (^1,^2; ^3,^4) 



dx\ dx 2 

d d 

-2(x 3 - x 4 ){- h — )£l_ 1 (x 1 ,x 2 ;x 3 ,X4) 



)[^(, 1 ,x 2 ;x3,x 4 )-(x3 + x 4 )^(x 1 ,x 2 ;x 3 ,x 4 ) 



-dxi dx 2 - 

F 4 (xt, x 2 , x 3 , x 4 ) = 2(lnx 4 - l)g (x 1 ,x 2 ) - 



6(x 3 — x 4 ) 6(0:3 In £3 — x 4 lnx 4 ) 



X\X 2 

x±x 2 - 2(xi + x 2 )(x 3 - x 4 ) 



X\X 2 



X\X 2 



, x 3 + x 4 dip 
fo{x 3 ,x 4 ) H — (x 3 ,x 4 ) 



x 4 (x 3 - x 4 ) gVo 

2 1^3) 2-4 J 



<9xj 



xix 2 dx 4 
+ x 4 ^2-)fi (xi,x 2 ; x 3 ,x 4 ) 



d d 2 
+ (- 1 + (x 3 + x 4 )-^ x 4 (x 3 - ^4)^3)^-1(^1, x 2 ; x 3 , x 4 ) 



d d \2 



dx 2 

f d d 



-dxi dx 2 



F 5 (xi,x 2 ,x 3 ,x 4 ) 



'<9xi (9X2 

-2(2 + lnx 4 )£ 01 (xi,x 2 ) 
x 3 - x 4 d<p 



fi 1 (xi,x 2 ;x 3 ,x 4 ) - (x 3 - x 4 )n o (xi,x 2 ;x 3 ,x 4 ) 

(XI 

dx. 

fi_ 1 (xi,x 2 ;x 3 ,x 4 ) 
1 



fi (xi, x 2 ; x 3 , x 4 ) - 2x 4 — — (xi, x 2 ; x 3 , x 4 ) 

ax 4 

2(x 3 - x 4 )y- h 



■(x 3 ,x 4 ) 



XiX 2 <9x 4 <9x 4 
d_ 
' dx 4 



-(p (x 3 ,x 4 ) 

X\X 2 



+ (l - (x 3 - x 4 )— )fi_i(xi,x 2 ;x 3 ,x 4 ) 

1 



F 6 (xi,x 2 ,x 3 ,x 4 ) = 2(2 + lnx 4 )^ 01 (xi,x 2 ) 

x 3 - x 4 dipo 



-ip Q (x 3 ,x 4 ) 



xix 2 

X1X2 ox 4 ox 4 
d 

1 - (x 3 - x 4 )— )fi_i(xi,x 2 ;x 3 ,x 4 ) . 



(Al) 



The concrete expression of $(x, y, z) can be found in ly, \2J\. In the limit z <C x,y, we 
can expand &(x,y,z) according z as 

z 2 z 3 z 4 

$(x, z) = <^ (x, y) + z<pi(x, y) + -y V2(x, y) + ^sOr, y) + ^a{x : y) 
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+2z(]nz-i)(l + Q ltl (x,y 

2 Anz 3w x + y 2xy yx 

Z V 9.\ A^(x-y) 2 (x-y) 3 ^ x> 



2z 3 /In 2 1 1 \ /_ 12xy 6xy(x + y) y- 



zz /mz LL \(-\ lzx 
-n\2\~W ~ 3fiA + (x. - 



+ 



(x — y) 2 ^ 3! 36^ ^ (x — y) 2 (x — y) 3 .< 
lf \nz 25 w 2x 3 + 58x 2 y + 58xy 2 + 2y 3 



-2z 

V 4! 288 A (x-y) 6 
24xy(x 2 + 3xy + y 2 ) ^ y x 
(x — y) 7 



^>o(z,y) = < 



(x + y)lnxlny + (x-y)0(x,y) , x > y ; 
2x In 2 x , x = y ; 

(x + y) lnx lny + (y — x)©(y, x) , x < y . 



<Pi{x,y) 



- In x In y - |±J 6(x, y) , x > y ; 
4 — 2 In x — In 2 x , x = y ; 

— lnx lny — |^6(y, x) , x<y. 



(2x 2 +6x;y) In x-(6xy+2;/ 2 ) In y _ 4xy q/ \ . 

(x— j/) 3 (x— tA 3 v ' "/ ' " ' 



V^2(x,y) = 



(a;-?/) 3 



2 — I — £_ ]ti T 

9x ^ 3x 111 X ' 

„2 



x = y ; 

(2x 2 +6x;y) In x-(6x;y+2;/ 2 ) lny _ Axy q/ \ , 

(x-j/) 3 r„_^3^(,y,x; , x y . 



(j/-x) 3 



12xy(x+y) 
(x-j/) 5 



0(x,y) 



2(x 2 +xj/+j/ 2 ) 



, 2(x 3 +14x 2 j/+llxi/ 2 ) lnx-2(y 3 +14xy 2 +Hx 2 j/) lny . 



(x-j/) 5 
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^2 lnx , 



150x 2 ^ 5x 2 
12xy(x+y) r\/ \ _ 2(x 2 +xy+y 2 ) 
(»-*) B U W> X ) (x-yY 



x = y; 



, 2(x 3 +14x 2 j/+llxj/ 2 ) lnx-2(y 3 +14xj/ 2 +llx 2 i/) lny 



(x-y) 5 



+ 



,y) 



' 48xy(x 2 +3xy+;/ 2 ) p ) / \ _ 2(3x 3 +61x 2 y+61xy 2 +3i/ 3 ) 

(X-J/) 7 yJ \- L - I U) ( X -yf 

4(x 4 +3x 3 y— 45x 2 y 2 -25xt/ 3 ) In x-4(t/ 4 +3y 3 x— 45x 2 y 2 -25j/x 3 ) lny 

(x-j/) 7 ' X > y 

2205x 3 + 210x 3 ^ n ^ ' x = Vi 

A8xy(x 2 +3xy+y 2 ) p.f \ _ 2(3x 3 +61x 2 y+61xy 2 +3j/ 3 ) 
(x-y) 7 (a;_ y )6 

4(x 4 +3x 3 j/-45x 2 j/ 2 -25xj/ 3 ) lnx-4(j/ 4 +3y 3 x-45x 2 t/ 2 -25j/x 3 ) In?/ 

' x < ^ • 
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Here, the function Q(x,y) is defined as 



2 

Q(x,y) =lnxln^-21n(x-y)ln^-2L« 2 (^) + — . (A8) 
x x x 3 
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